Theory of electromagnon in the multiferroic Mn perovskites 
— Vital role of higher harmonic components of the spiral spin order — 
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We study theoretically the electromagnon and its optical spectrum (OS) of the terahertz-frequency 
regime in the magnetic-spiral-induced multiferroic phases of the rare-earth (_R) Mn perovskites, 
7?Mn03, taking into account the elliptical deformation or the higher harmonics of the spiral spin 
configuration, which has been missed so far. A realistic spin Hamiltonian, which gives phase dia- 
grams in agreement with experiments, resolves a long standing puzzle, i.e., the double-peak structure 
of the OS with a larger low-energy peak originating from magnon modes hybridized with the zone- 
edge state. We also predict the magnon branches associated with the electromagnon, which can be 
tested by neutron-scattering experiment. 

PACS numbers: 75.80.-|-q, 75.40.Gb, 75.30.Ds, 76.50.+g 



The charge dynamics below the Mott gap in Mott in- 
sulators is an issue of intensive recent interest. The rich 
structures of the low-energy optical spectrum (OS) are 
associated with the spin degree of freedom in so-called 
multiferroics, which shows both magnetic order and fer- 
roelectricity The spontaneous electric polarization 
in these materials is driven by the magnetic ordering, 
and the strong coupling between the dynamics of electric 
polarization and magnetism is inevitable. In particular, 
in the rare-earth (R) perovskite manganites i?Mn03 0, 
which is the main target of this paper, the relativistic 
spin-orbit interaction and the spin current play essen- 
tial roles in the multiferroic behavior Excitations 
of the magnet are usually described as the spin wave or 
the magnon, i.e., a harmonic oscillation of spins around 
their ground-state configuration. It couples to the elec- 
tric polarization in the multiferroics, and thus is termed 
electromagnon @. It is crucial to understand the OS 
of the electromagnon to design the giant magnetoelectric 
coupling in the terahertz-frequency regime. 

In i?Mn03, frustration between the spin-exchange in- 
teractions leads to a non-collinear spin spiral in the 
ground state. This simple idea, however, cannot explain 
the rich phase diagrams in the plane of Mn-O-Mn bond 
angle and temperature, and also those under magnetic 
fields. In real systems, there arc other interactions origi- 
nating from the relativistic spin-orbit interaction such as 
the magnetic anisotropy and the Dzyaloshinskii-Moriya 
(DM) interaction. By studying a realistic spin Hamil- 
tonian taking into account these interactions, the phase 
diagrams including the spin-flop transition have been un- 
derstood except for the collinear E-type spin phase 0- 
It turned out that coupling of the spins to phonons given 
by Hsp = X^ij gQijSi ■ Sj plays an important role to sta- 
bilize the E-type phase, which is expected to play some 



roles even in the neighboring spiral phases. 

On the other hand, after the first experimental obser- 
vation of the terahertz OS in the multiferroic -RMnOs [1, 
Q, it was interpreted as a collective mode correspond- 
ing to rotation of the spin-spiral plane associated with 
fluctuations of the polarization direction [lo| . Later it 
turned out that the selection rule and the magnitude of 
oscillator strength ruled out this interpretation [ll|, [l^] , 
and a new mechanism has been searched for. The most 
promising candidate is the conventional magnetostriction 
mechanism 



given by 



13l . Il4j , where the electric polarization P is 
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Here the vector Ily is nonzero in i^MnOa since the in- 
version symmetry is absent at the center of the Mn-O- 
Mn bond because of the orthorhombic lattice distortion 
and/or the staggered d^^2_r2 /d^y2_r2 orbital ordering. 
This contribution cancels out in the ground state due to 
the symmetry, but the dynamical fluctuations of P con- 
tribute to the electromagnon excitation. Especially, in 
the non-collinear ground state, the single magnon pro- 
cesses at the zone edge originate from Eq. ([T]). However, 
this scenario cannot explain the low-energy peak at ^2- 
3 meV (see insets of Fig. [5]), which is comparable to or 
even larger than the high-energy peak at ~5-8 meV in 
DyMnOg [llj, TbMnOg [H, and Eui^^Y^MnOg 0. 
Therefore, the puzzle still remains. 

A clue to this issue is the proximity to collinear spin 
phases, i.e. the A- type and E-type spin phases. Near the 
phase boundary, the spin configuration is not a simple 
spiral but is subject to the significant elliptical modula- 
tion and contains higher harmonics [l3|, which is sensi- 
tively enhanced by the tiny spin-phonon coupling or by 
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FIG. 1: (Color online) (a) Superexchange interactions in 
_RMn03 described by the Hamiltonian Eq. (2). (b) Mod- 
ulations of the in-plane nearest-neighbor ferromagnetic ex- 
changes under E\\a. (c) Those under E\\b. Here FM and 
AFM denote ferromagnetic and antiferromagnetic exchanges, 
respectively. 



the weak magnetic anisotropy. In this paper, we study 
the role of this higher harmonics on the electromagnon 
excitation, and resolve the long standing puzzle of the 
OS in the terahertz-frequency regime. 

We start with a realistic spin model for the Mn per- 
ovskites • The model is basically a classical Heisenberg 
model on a cubic lattice, which contains the frustrating 
spin exchanges, the DM interaction, the single-ion spin 
anisotropies, and the biquadratic interaction. The Mn 
spins are treated as classical vectors. The Hamiltonian 
consists of five terms as 7^ = 'Hex + + Tif^^ + Hdm + 
■Hbiq, with 
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i-th Mn 



where 

ion with respect to the cubic x, y and z axes. For the a, 
b and c axes, we adopt the Pbnm notation [see Fig. [TJa)]. 
The first term Hex describes the superexchange interac- 
tions as shown in Fig. [TJa). The frustration between fer- 
romagnetic Jab and antiferromagnetic J2 in the ab plane 
results in a spiral spin order, while the inter-plane anti- 
ferromagnetic Jc causes a staggered stacking along the c 
axis. The terms Hf? and stand for the single-ion 

old bid ' 



anisotropies. "H^g^ makes the magnetization along the c 
axis hard. On the other hand, Ti^^^ causes an alternate ar- 
rangement of the local hard and easy magnetization axes 
in the ab plane due to the staggered orbitals. The term 
Hdm denotes the DM interaction. The DM vectors dij 
associated with Mn-O-Mn bonds are expressed using five 
DM parameters as given in Ref. because of the crystal 
symmetry; aab, Pab and ^ab for dij on the in-plane bonds, 
while ac and j3c for dij on the inter-plane bonds. The last 
term T^biq represents the biquadratic interaction work- 
ing between adjacent two spins in the a6-plane, which 
originates from the spin-phonon coupling [19[. We have 
microscopically determined the parameter values 0]. 

We perform calculations using two sets of the model 
parameters (A and B) as (A) Jab=:— 0.74, J2=0.74, 
Je=1.2, (aa&,/3a&,7a&)=(0.1, 0.1, 0.16), (a„/3,)=(0.4, 
0.1), L>=0.24, E=Q.i, and Bbiq=0.025, and (B) 

Ja6 = -0.7, J2=0.96, Jc = 1.0, (aafcAfc,7a6) = (0.1, 0.1, 

0.12), (ac,/3c)=(0.45, 0.1), £1=0.2, E=Q.2b, and 
^biq=0.025. Here the energy unit is meV. These pa- 
rameter sets give the a6-plane spin spiral propagating 
along the b axis with six times periodicity (qi,=Tr/3) and 
the 6c-plane one with five times periodicity ((7b=27r/5), 
respectively. The former spin state resembles the ab- 
plane spiral in Eui_j;Yj;Mn03 (x=0.45) with qh^O.Sn, 
while the latter resembles the 6c-plane spiral in DyMnOa 
with (7t,=0.397r. Note that we adopt the commensurate 
spin states for convenience of the finite-size calculations, 
whereas the actual spin states are incommensurate. How- 
ever, conclusions of this paper are never affected by this 
difference. Sizes of the systems used are 20x20x6 and 
18x18x6 for respective cases, which match the periodic- 
ities of spirals. 

We study the electromagnon and magnon excitations 
by numerically solving the Landau-Lifshitz-Gilbcrt equa- 
tion using the fourth-order Rungc-Kutta method. The 
equation is given by 



dt 



H 



off 



s 



where q;g(=0-1-0.2) is the dimensionless Gilbert- 
damping coefficient. We derive an effective local mag- 
netic field H![^ acting on the i-th Mn spin Si from the 
Hamiltonian H as H!?^ = —dH/dSi. Considering the 
strong reduction of the Mn magnetic moment revealed 
by neutron-scattering experiment [l7j . we set |5| = 1.4. 

For the origin of the electromagnon excitation, we con- 
sider the coupling ■ P between the external electric 
field E and the spin-dependent electric polarizations P 
given by Eq. ([1]). Noticeably this coupling effectively 
modulates the nearest-neighbor ferromagnetic exchanges 
in the ab plane from JabSi ■ Sj to {Jai, + E ■ Tlij)Si ■ Sj. 
More concretely, the application of E\\a [E\\b] corre- 
sponds to modulations of the in-plane spin exchanges 
shown in Fig. (Ub) [Fig. (He)]. 

We apply the electric field E\\a oi E\\b a.s a short pulse 
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FIG. 2: (Color online) Calculated electromagnon optical spec- 
tra for (a) a6-plane spiral state (gi,=7r/3) with parameter set 

A, and (b) foc-plane spiral state (qf,=27r/5) with parameter set 

B. Values of qg used for the calculations are 0.1 and 0.2 for (a) 
and (b), respectively. Insets show the experimental spectra 
of Euo.55Yo.45Mn03 [iil and DyMnOs [Tj, respectively. 



at t=0. Then we trace the time evolution of P given 
by Eq. (1). The electromagnon spectrum, Im e(a;), is 
calculated from the Fourier transformation of P{t). 

In Fig. [21 we depict calculated electromagnon spectra 
for the cases of (a) 6c-plane spiral state with parameter 
set A and (b) a&-plane spiral state with parameter set B. 
Irrespective of the spiral-plane orientation, large spectral 
weight emerges at low energy when E\\a. In contrast, we 
observe no response to E\\b for both cases in agreement 
with the experiment [12[ . For comparison, we display the 
experimental spectrum of Euq. 55 Yo.45Mn03 [l^ and that 
of DyMnOs [11| in the insets. For both cases, we obtain 
fairly good agreement between theory and experiment. 
In the following, we discuss the results for the a6-plane 
spiral case. Similar discussion can be repeated for the 
&c-plane spiral case. 

To identify origin of the two-peak structure of the elec- 
tromagnon spectrum, we calculate electromagnon spec- 
tra and magnon-dispersion spectra along k = (0,fcb,0), 
for various cases of interactions (see Fig.[3l). The magnon 
spectra are calculated from the Fourier transformation of 
the space- and time-domain simulation data for the spin 
dynamics of (5511 (a&-spiral plane) after applying H as a 
short pulse to a single-site spin in the a6-spiral ground 
state. When the Hamiltonian consists of only the spin- 
exchange term and the single- ion anisotropy D term, i.e., 
Hex and "Hja [see Fig. Elja)], the spin spiral has a uni- 
form rotation angle. In practice, for the case of pure 
spiral order, translational symmetry is conserved upon 
the one lattice-unit translation if it is accompanied by an 
appropriate rotation of the spin axes. Then despite the 
long-period magnetic structure, matrix elements which 
mix magnon branches do not exist. In this case, we can 
see only one peak at a rather high energy of ~8 meV 
corresponding to that in Fig. |31Ja). Incorporation of the 
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FIG. 3: (Color online) Calculated intensity map of magnon 
(left panel) and electromagnon OS (right panel) for each 
Hamiltonian with successively adding the interactions, (a) 
'Hcx+'HSu giving a spin spiral with a uniform rotation angle, 
(b) Adding "Hdm induces negligibly small changes, (c) Incor- 
poration of Tif^^ causes the elliptical deformation or higher 
harmonics of the spin spiral, resulting in the magnon foldings 
and evolution of the lower-energy peak in the OS. (d) For full 
Hamiltonian including also "Hbiq, the lower-energy peak in the 
OS is further enhanced. 



DM interaction T^dm does not change the spectral shape 
[see Fig.Hb)]. 

Further adding the single-ion anisotropy E term 
gives rise to folding and anticrossing of the magnon dis- 
persions. Namely, in the extended Brillouin zone picture, 
magnon branches separated by the reciprocal lattice vec- 
tor G=(0, Qb, 0) are mixed with each other where qi,=Tr/S 
in the present case. This gives rise to changes in the spec- 
tral shape, i.e. another peak of the OS at a lower energy 
of ^3 meV appears, and the higher-lying peak slightly 
splits into two peaks as shown in Fig. ^c). Finally we 
can sec dramatic enhancement of the lower-lying peak 
in the OS when the Hamiltonian is full including the bi- 
quadratic term Hbiq as shown in Fig. El^d). Comparing 
the electromagnon spectra with the magnon spectra at 
the zone edge, i.e. fe=(0, tt, 0), we indeed confirm that 
the electromagnon corresponds to the magnon modes hy- 
bridized with the zone-edge state due to the form factor 
of the magnetoelectric coupling in Eq. ([1]) as described 



in Ref. jl3|. Strong magnon anticrossing causes rather 



fiat magnon dispersions. The magnon branches predicted 
here would be observed experimentally. 

We can see vital roles of the single-ion anisotropy 
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FIG. 4: (Color online) Calculated electromagnon spectra for 
various strength of (a) single-ion anisotropy E and (b) bi- 
quadratic interaction -Bbiq- Insets show calculated ellipticity 
of the afe-plane spin spiral (ma/mi,) as functions of E and 
-Bbiq, respectively. 



and the biquadratic interaction Ubiq in Figs, mja) and 
m^b), which show the calculated spectra for various val- 
ues of the coupling constants E and -Bbiq, respectively. 
Surprisingly the lower-lying peak is enhanced strongly by 
the weak anisotropy or by the tiny biquadratic interac- 
tion. In the insets of Figs.[3ja) and El^b), we show calcu- 
lated ellipticity of the spin spiral for each case, which is 
defined as the ratio of amplitudes between the a-axis spin 
component and the 6-axis spin component [ma/rrib) |20|. 
We see that the enhancement of the low-lying peak are 
accompanied by the decrease of ellipticity. Note that the 
elliptical deformation of the spin spiral is more significant 
when the ellipticity is decreased more from unity. This 
indicates that the origin of the low-energy peak at ~2-3 
mcV is the magnon-branch mixing due to the elliptical 
modulation of the spin spiral. 

Let us now discuss the R dependence of the electro- 
magnon spectra of /JMnOa. In TbMnOs, the spectrum 
has more weight at the higher-energy peak |15| . while 
in DyMnOs with a smaller ionic i?-site radius {rp), the 
spectral weight is shifted to lower energy 
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In the 

systematic study of Eui_.^Y2.Mn03, transfer of the spec- 
tral weight to lower energies has been reported as the 
Y content x is increased (or as the averaged is de- 
creased) [31 . Note that with decreasing r^, we approach 
the E-type phase, which results in a stronger influence of 
the biquadratic interaction. So far, the magnon spec- 
trum in i?Mn03 with a sp iral spin order has only been 
reported for TbMnOs [2l|. In the compounds closer to 
the E-type phase, we expect that the magnon spectrum 
exhibits significant foldings. Further experimental data 
for these compounds are necessary to examine our theory. 

In some of the i?Mn03 compounds, deviation of fre- 
quencies between the magnetic resonance and the electro- 
magnon has been reported [1^. Within our theory, the 
frequency of the magnetic resonance at the zone center 



fc=0 coincides with that of the electromagnon at fc=(0, 
TT, 0) only if they are connected with a multiple of G=(0, 
qb, 0). When the spiral magnetic order is incommensu- 
rate, magnetic resonance and electromagnon absorption 
occur at different points in the Brillouin zone. Thus fre- 
quencies of these two spectra do not necessarily coincide 
with each other. 

In summary, we have studied the electromagnon dy- 
namics taking into account all the relevant interactions 
and anisotropics in the spin Hamiltonian for KMmO^. It 
is found that the elliptical modulation or higher harmon- 
ics of the ground-state spiral spin configuration has cru- 
cial influences on the foldings of the magnon dispersion 
and also the electromagnon spectrum, which resolves the 
long standing puzzle of the low energy peak in the OS. 
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